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摘 要 : 为 了 改善 粒子 群 优化 算法 的 求解 性 能 ， 提 出 了 一 种 基于 单纯 形 搜索 和 粒子 群 优化 的 混合 算法 。 该 算法 一 方面 
自 适应 地 确定 惯性 权重 、 认 知 以 及 社会 参数 来 达到 免 参 数目 的 , 另 一 方面 利用 单纯 形 搜索 来 引导 部 分 粒子 的 搜索 方向 ， 
从 而 加 速算 法 收敛 。 数 值 实验 结果 表明 ， 与 传统 的 粒子 群 算 法 和 其 他 基于 单纯 形 的 粒子 群 算法 相对 比 ， 该 算法 在 评估 
次 数 、 求 解 精度 方面 表现 良好 。 

关键 词 : 直接 搜索 ; 单纯 形 搜索 ; 粒子 群 优化 

中 图 分 类 号 : TP301.6 doi: 10.3969Vj.issn.1001-3695.2018.06.0465 


Particle swarm optimization based on simplex search 
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Abstract: In order to improve the performance of particle swarm optimization, the paper proposed a hybrid algorithm based on 
simplex search and particle swarm optimization. On the one hand, the algorithm determines the inertia weight adaptively, 
cognition and social parameters for the purpose of avoiding parameters. On the other hand, the simplex search is used to guide 
the direction of several particles, thus the convergence of the algorithm is accelerated. The results of numerical experiments 
show that the proposed algorithm has better performance than the compared algorithms in terms of function evaluations and 
accuracy. 
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果 ， 并 在 工程 领域 到 广泛 应 用 [5~!。2007 年 ，Fan 和 Zahira 提 
出 了 一 种 求解 无 约束 优化 问题 的 混合 算法 ， 该 算法 将 单纯 形 搜 
粒子 群 算法 是 由 Kennedy 和 Eberhart 于 1995 年 提出 ,， 它 。 ” 索 与 粒子 群 算法 相 结 合 ， 以 提高 粒子 群 优 化 的 收敛 速率 ， 其 采 
通过 模拟 鸟 群 郧 食 过 程 中 的 迁徙 和 群 聚 行为 来 寻找 全 局 最 优 解 有 的 算法 结构 为 后 续 算法 设计 提供 了 思路 ， 但 该 算法 参数 设置 
[3] 。1998 年 Shi 等 人 引入 了 惯性 权重 机 制 ， 使 得 粒子 群 算法 。 较 多 ， 且 对 较 高 维 数 的 优化 测试 函数 的 求解 性 能 有 待 进一步 提 
能 得 到 进一步 改善 m 。 与 传统 优化 方法 相 比 ， 粒 子 群 算 ”高 09。2008 年 ，Hsu 和 Gao 在 粒子 群 优 化 中 嵌入 了 单纯 形 搜索 
法 具有 如 下 优点 : 不 需要 目标 函数 导数 等 解析 信息 ， 有 具有 很 好 机制 ， 该 机 制 每 隔 上 次 粒子 群 优化 算法 迭代 便 执 行 1 次 单纯 形 
I 通用 性 ; 不 易 陷 入 局 部 最 优 , 具有 较 强 的 鲁 棒 性 ; 流程 简单 ， 搜索 ， 从 而 不 但 可 以 减少 函数 评估 次 数 ， 而 且 较 好 地 平衡 了 搜 
易于 实现 。 然 而 ， 粒 子 群 算法 的 启发 式 机 制 往往 使 得 算法 易于  ” 索 和 探索 之 间 的 矛盾 ， 但 是 数值 实验 部 分 几乎 只 评估 了 2-5 维 
早熟 ,而 且 对 于 特定 问题 ,算法 性 能 过 多 依赖 于 对 参数 的 设 定 。 化 函数 的 求解 性 能 071。2010 年 ， 任 小 康 等 人 提出 了 一 种 基于 
如 今 许 多 改进 算法 也 未 能 很 好 解决 参数 设置 过 多 ,， 且 参数 设置 ”单纯 形 法 的 量子 粒子 群 优 化 算法 ， 该 算法 先 用 量子 粒子 群 算法 
依靠 经 验 &%I0， 而 对 于 简洁 的 无 参数 改进 算法 而 言 ， 算 法 易 早 “进行 全 局 搜索 ， 然 后 依据 种 群 适应 度 标准 差 确定 是 否 利用 单纯 
熟 收 么 0 。 单 纯 形 搜索 能 够 利用 问题 本 身 较 弱 的 解析 性 质 ， 形 法 进行 局 部 搜索 Q81， 但 是 如 何 根据 种 群 适应 度 标准 差 判断 种 
理论 基础 丰富 ， 搜 索 具 有 针对 性 ， 收 敛 速 度 较 快 ， 并 且 和 粒子 半 是 否 陷 入 “早熟 ”是 个 困难 的 事情 ， 从 而 限制 了 算法 的 进 
群 优化 有 着 天 然 的 共性 ， 便 于 取长补短 。 考 虑 到 这 两 个 算法 的 。” ” 步 推广 。 2018 年 , 武 可 栋 等 人 在 粒子 群 优化 算法 的 初始 化 过 程 
特点 ， 出 现 了 许多 基于 单纯 形 的 粒子 群 算法 21 。 中 加 入 了 单纯 形 搜索 ， 利 用 单纯 形 法 对 初始 值 进 行 预 处 理 09， 
近年 来 ， 有 许多 基于 单纯 形 搜索 的 粒子 群 优 化 算法 研究 成 ”并 针对 20 个 10 维 测试 函数 进行 了 数值 模拟 ， 实 验 结果 表明 其 
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相 比 标准 粒子 群 优化 算法 具有 更 高 的 求解 精度 ， 
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但 是 该 算法 存 


在 参数 设置 困难 ， 评 估 次 数 较 多 的 缺点 。 

为 了 改善 粒子 群 优化 算法 的 求解 性 能 ， 本 文 提出 了 一 种 基 
于 单纯 形 搜索 的 粒子 群 优化 算法 ， 该 方法 一 方面 利用 类 电磁 机 
制 原理 023 自 适应 地 确定 惯性 权重 、 认 知 以 及 社会 参数 来 达到 免 
参 目的 , 另 一 方面 利用 单纯 形 搜索 来 干预 每 次 搜索 方向 及 步 长 ， 
从 而 加 速算 法 收敛 。 


1 “粒子 群 算 ; 


粒子 群 算法 (particle swarm optimization，PSO) 的 思想 来 
源 于 对 鸟 群 飞行 或 驶 食 过 程 中 的 迁徙 和 群集 行为 0 7339。 鸟 仅 
仅 是 追踪 它 有 限 数 量 的 邻居 ， 但 最 终 的 整体 结果 是 整个 鸟 群 好 
像 在 一 个 中 心 的 控制 之 下 ， 即 复杂 的 全 局 行为 是 由 简单 规则 的 
相互 作用 引起 的 。 不 同 于 达尔 文 “ 适 者 生存 ， 优 胜 劣 汰 ”进化 
思想 ， 粒 子 群 优化 算法 是 通过 个 体 之 间 的 协作 来 寻找 最 优 解 。 

PSO 算法 描述 如 下 : 

a) 初 始 化 种 群 P(1) ,速度 V(7) , 认 知 系数 co, 社 会 系数 c2， 
惯性 权重 w， 令 410; 

b) 计 算 每 一 个 粒子 X, e P(1) 的 适应 值 , 个 体 经 历 最 优 位 置 
g,(t)， 种 群 全 局 最 优 位 置 G(7); 

中 更 新 每 一 个 粒子 X,eP(1) 的 速度 Vj,(1+1) 和 位 置 
X,(t+1): 

for X,(t)eP(?) 


for j=1 ton 
V(t+D)=w:V,() 
+ci:rand (8 ,(D) —X, iD) (1) 
+c :rand:(G,(D) —X,,()) 
end 
X(t+D)= XD +V+D) (2) 
end 
gd) 若 满足 终止 条 件 ， 停 机 否则， =r+1， 转 b)。 
这 里 ，V, (1) 表示 VD) 的 第 j 个 分 量 ，X (7) 表示 X,() 的 
第 j 个 分 量 ，n 表示 维 数 ， 即 分 量 个 数 ， cl 和 cz 是 两 个 正常 
数 ，w 是 惯性 权重 ，rand 是 一 个 (0,1) 之 间 的 随机 数 。 
在 步骤 c) 中 , 式 (1) 说 明 粒 子 的 速度 取决 于 上 次 的 速度 更 
新 ， 它 自己 所 到 达 的 历史 最 佳 位 置 和 全 种 群 中 所 有 粒子 的 历史 
最 佳 位 置 ; 式 (2) 表 明 每 个 粒子 的 位 置 在 搜索 域 中 的 更 新 方式 。 


2 单纯 形 搜 索 


单纯 形 (simplex) 是 代数 拓扑 中 的 基本 概念 ， 是 三 角形 和 
四 面体 的 一 种 泛 化 。n 维 单纯 形 是 一 个 由 n+1l 个 点 组 成 的 凸 包 。 
例如 : 1 维 单纯 形 就 是 线段 ，2 维 单纯 形 就 是 三 角形 ; 3 维 单纯 
形 就 是 一 个 四 面体 。 

单纯 形 搜 索 (Nelder-mead simplex search，,， 简称 NM) P29 通 
过 反射 、 扩 展 、 压 缩 等 一 系列 初等 几何 变换 操作 来 寻求 问题 的 
最 优 解 。 每 次 变换 后 ， 试 图 用 一 个 更 好 的 顶点 蔡 换 当前 已 知 最 


Im 


差点 。 
以 最 优化 问题 min ,了 (X) 为 例 , 分 别 计算 3 个 顶点 Xi ， 
2 > Xi 的 隔 数值 fiign 2? foe ? fi 用 不 妨 假 设 ， fiign > Tg > fy 


基本 操作 (图 1~4) 如 下 : 
a) 反射 。 取 久 


和 X 的 中 点 xX， 计算 反射 点 


low cent 


Xs = bE + (Xo _ Xion ) 
b) 扩展 。 著 fo。>f。， 则 计算 扩展 点 
有 二 Xo + 2 (X,, A ) 


c) 压缩 。(@) 著 fi >fs>f.。， 计 算 向 外 压缩 点 


Xo = Xoo +0.5*(X,, — Xe) 


cont cent 


， 用 多 ,取代 义 


于 
石 fro 2 cont cont high » 


(b) 若 fj>fim ， 计算 向 内 压缩 点 


xX — +0.5* (Xo —X.,) 


cont cent 


br 取代 X pion 由 


一 fhigh ; 则 向 最 好 点 和 
BW 对 Vi, i low, X=05*(X,+ Xi,) : 


名 fhigh > fa S 
d) 收 缩 若 / 


cont 


收缩 其 余 顶 点 ， 


low 


Xsec Xo 
图 1 反射 操作 图 2 扩展 操作 
Fig.l Reflection Fig.2 Expansion 
Xion Xiow 
Xigh Xjow 


Xeo a 
X 
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(a) 向 外 压缩 操作 


(a) Outside contraction 


(b) 向 里 压缩 操作 
(b) Inside contraction 
图 3 压缩 操作 


Fig.3 Contraction 
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图 4 收缩 操作 
Fig.4 Shrink 
NM 算法 描述 如 下 : 


a) 计 算 n 维 单纯 形 的 ntl 顶点 X 的 函数 值 f ,i=1…,n 
选取 最 差 解 Xiw ， 次 差 解 X,。 ， 最 好 解 X,,， 中 点 


XX =0.3*(Xw+X) ， 及 其 对 应 函数 值 为 fi ,f， fi,， 
et 9 
b) 计 算 反 射 点 X=X +(Xw 一 Xi ) 以 及 对 应 函数 值 
为 f.,。 
9 进行 单纯 形变 换 : 
if ju > fo 
if fr = 3 ss 
计算 压缩 点 Xn = > 二 0.5 和 (Xiion — X oi) 
如 果 fe < fis ” 则 Xion Sr ; 口 否则 向 Xj 进行 收 


缩 操作 。 
elseif fi > fry > fee 


和 =X +0.5* (X | 


cont cent 


如 果 六 < fw ， 则 Xi = XX; 否则 向 X ,进行 收缩 
操作 。 
end 
> 
else 
计算 扩展 点 X ,=X +2*(X,y 一 Xo) ， 如 果 f, < 


Pa ， 则 Xjion =X。, ; 否则 ， Xigh = Xo 3 


end 

通过 步骤 c) 的 操作 ,新 单纯 形 的 %41 个 顶点 中 至 少 存在 一 
个 顶点 进行 了 改变 。 
3 ”基于 单纯 形 搜索 的 粒子 群 优化 算法 
3.1 PSO 算法 改进 
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参考 点 X, (1) ，X,,(?) ， 基 点 X.;(D 为 任意 选取 的 其 他 三 个 
不 同 的 粒子 ， 即 irlr2zr3.。 

式 (5) 采用 了 贪 禁 选 择机 制 更 新 粒子 X,(t) ， 即 ， 仅 仅 当 
候选 粒子 TCD 的 评估 值 f(T() 小 于 粒子 X (CD 的 评估 值 f(7) 
时 ， 用 候选 粒子 蔡 换 掉 xX,(1)。 式 (4) 是 利用 基点 X.,(D 以 及 
其 与 参考 点 之 间 的 信息 共享 来 构造 候选 粒子 T(1) 。 

记录 了 参考 点 X,() 与 基点 XX,,(n) 之 间 的 交流 信息 。 
f(D < f(D 时 ,5 2[FGX OO)-FXOD) 71D>0， 古 将 引导 
XX,3(7) 向 XGO) 方向 搜索 ， 而 且 ， 评 估 值 差异 越 大 ， 在 方向 
(X,(D) 一 X,s(D)) 上 搜索 步 长 8, 就 越 多 , 反之 , 差异 越 小 ，5 就 
越 多 ; 当 fj(1) > f(D) 时 ，S,。 <0， 厂 将 引导 XX,,(7?) 在 ae 
内 搜索 ， 这 时 搜索 步 长 9, 刻画 了 邻 域 的 大 小 ,而且 ， 评 估 值 差 
异 越 大 , 邻 域 9, 就 越 小 , 反之 , 差异 越 小 , 邻 域 就 越 大 。 因此 ， 
通过 这 种 方式 ， 一 方面 可 以 很 好 地 平衡 局 部 搜索 和 全 局 探索 之 
间 的 矛盾 ， 另 一 方面 算法 自 适 应 地 在 加 速 收敛 和 提高 种 群 多 样 
性 之 间 转 换 。 
司 样 的 处 理 思 路 ，F1 记录 了 参考 点 X,,() 的 个 体 经 历 最 优 
gw(1) 与 基点 X,,(?) 之 间 的 信息 。 
姑 为 必 有 f,(1) < 请 (D ， 所 以 到 引导 基点 X,(D 向 种 群 
全 局 最 优 CD 方向 搜索 。 

这 种 引入 两 个 粒子 间 的 信息 交流 机 制 类 似 于 电磁 机 制 。 把 
每 个 粒子 都 当 作 一 个 点 电荷 ， 函 数值 越 好 《小 ) 点 电荷 电量 越 
高 ， 对 其 他 粒子 的 吸引 力 越 强 。 一 个 粒子 的 下 一 步 搜索 方向 是 
其 他 粒子 对 它 施 加 的 “合力 ”， 因 此 这 种 机 制 不 仅仅 有 利于 改善 
算法 性 能 0220， 而 且 能 够 达到 免 参数 目的 。 

3.2 算法 描述 

a) 初 始 化 种 群 ， 种 群 规 模 POP=3n+1， 变 异 系数 Cr=0.1， 
t=0; 

b) 排 序 。 
序 ; 

ce) 单纯 形 搜索 。 将 前 n+1l 个 最 好 的 粒子 作为 单纯 
执行 NM 算法 ， 改 进 n+1l 个 顶点 中 最 差 粒 子 ; 

qd) 粒 子 群 优化 。 以 剩余 的 2n 个 粒子 构成 粒子 种 群 PLD)， 对 


I 


计算 函数 值 f(z)， 并 按照 函数 值 从 小 到 大 进行 排 


NN 


顶点 ， 


在 传统 粒子 群 优化 中 ， 第 i 个 粒子 的 位 置 更 新 基于 三 个 医 
素 ， 见 式 (1) (2)。 然 而 ， 在 基于 单纯 形 搜 索 的 粒子 群 优 
化 算法 (Nelder-mead simplex particle swarm optimization，NM- 
PSO) 中 , 对 于 第 t 代 种 群 P(r) 中 第 i 个 粒子 X,(1) ,通过 式 (3) 
~ (5) 更 新 : 


V=R+E+b (3) 
TD)=X,D+V (4) 
X(t+D=T@ 这 f(TOD)<fD (5) 
这 里 : 
R=(X,(D)— XD) [FX —f(XD)/D 
Fi=(g,2(D)— XD) [F(XD)-f(8,D))]/D 
EE=(G@ -xX,D)*[f(X,D) -GDND 
D= fi — fo 


每 一 个 粒子 X,(1) e PCD) ， 
(a) 按 照 式 (4) 计算 候选 粒子 T(1)=X(D)+h+R+b; 
(b) 边 界 检查 。 对 vj ， 如 果 了 DOD<ZB， 则 TD=ZB: 如 
果 T (D>UVB8B,， 则 TT,(D=UB; 
(9] 贪 禁 选 择 。 如 果 f(T(D) < f(D， 则 X(t1+1)=T(7)， 
否则 X,(1+1)=X,(7); 
(qd) 变异 。 对 vj ， 如 果 rand> Cr， 则 XIC+D=X CD 
e 如 果 满 足 终止 条 件 , 则 停机 ; 否则 ，r =z+1， 转 步 又 b)。 
为 了 增加 种 群 的 多 样 性 ， 避 免 算法 过 早 陷 入 局 部 最 优 ， 在 
步骤 c 中 ，NM 算法 只 执行 一 步 迭 代 ， 步 又 四 中 ， 增 加 了 变异 
操作 。 


4 ”数值 实验 


胡 锦 帆 ， 等 : 


4.1 对 照 算法 介绍 

本 文 将 比较 包括 本 文 算法 在 内 的 五 个 算法 的 数值 实验 结果 。 
首先 介绍 其 他 四 种 对 照 算法 。 

1) 惯性 权重 递减 粒子 群 算法 (degressive inertia weight 
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随机 数 ， 即 w=rand (0.5,1) 器。 
3 ) 混合 单纯 形 粒子 群 算法 09 
optimization, HSPSO) 
HSPSO 算法 在 每 一 次 迭代 中 ， 


其 余 步 又 均 与 PSO 保持 一 


致 。 


hybrid simplex particle Swarm 


方面 利用 单纯 形 法 对 最 好 


的 前 N+1 个 粒子 进行 更 新 , 另 一 方 


面 利用 粒子 群 优 化 算法 对 剩 


余 的 2N 粒子 进行 改进 。 在 单纯 


NY 


法 每 一 次 执 


行 中 ， 只 需 和 迭代 


particle swarm optimization, DPSO) 


在 标准 


粒子 群 算法 中 ，{ 


次 速度 更 新 的 影响 力 相 比 另 
当 w 值 较 大 时 ， 


Y (0) 的 范 数 较 大 ， 每 次 


久 性 权重 w 对 WV , (71) ， 
外 两 个 分 量 


即 该 粒子 上 
5 据 主 要 位 置 ， 所 以 ， 


区 


和 迭代 粒子 的 “ 


离 ” 较 长 ， 有 利于 全 局 搜索 


T 


以 找到 最 优 解 所 在 的 邻 


当 w 较 小 时 ， 


力 回 。 因 
前 期 全 局 搜索 能 力 加 
能 够 使 算法 整体 性 能 提高 。 


三 个 分 量 对 迭代 方向 影 
相对 较 小 的 区 域内 能 够 更 精确 地 搜索 ， 即 加 强 了 
比 ， 如 果 使 惯性 权重 w 能 多 


=I 


力 相 近 ， 


强 , 在 运 


惯性 权重 w 的 递减 公式 如 下 : 


前 进 距 
或 。 反 之 ， 
时 算法 在 一 个 
局 部 搜索 的 能 
多 随 夫 代 次 数 递 减 ， 则 算法 


行 后 期 局 部 搜索 能 力 加 强 ， 


w= Wax | & LS ] * (Wax Wnin ) 
ZIE7 
其 中 : w。 ，w, 分别 为 设 定 的 w 的 最 大 值 和 最 小 值 ，Shi 等 人 
推荐 这 两 个 参数 的 值 分 别 为 1.4 和 01， iier ， ifer 分别 为 当 


el 


选取 前 nt1 好 的 粒子 
构成 单纯 形 ， 改 良 单 
纯 形 中 最 差 粒子 


图 5 NM-PSO 流程 图 


前 迭代 次 数 和 最 大 迭代 次 数 。 


初始 化 种 群 ， 随 机 
生成 3xn+1 个 粒子 


种 群 按 照 函数 
值 排序 


剩 下 较 差 的 2*n 个 粒 
子 进行 粒子 群 算法 改 


进 


Fig.5s Flowchart of NM-PSO 


采用 此 算法 作为 比较 算法 之 一 的 原因 是 ， 此 算法 也 采用 自 
适应 参数 的 机 制 ， 只 是 这 个 自 适应 机 制 较为 简单 ，NM-PSO 的 
自 适应 机 制 较为 复杂 ， 经 过 对 比 能 够 显现 出 算法 在 使 用 自 适 应 


参数 机 制 由 简单 到 精密 
2) 随机 惯性 
Swarm optimization，RPSO ) 


RPSO 是 将 PSO 中 的 惯 


Li 
T 


的 性 


能 提升 。 
权重 粒子 群 算法 (random inertia weight particle 


性 权 重 w 


改 为 [0.5,1] 中 的 均匀 分 


布 


次 ， 便 


新 的 粒子 蔡 换 原单 纯 


的 最 差 顶 点 ， 这 样 不 仅 可 以 


加 RS 


大 大 降低 算法 的 计算 复杂 度 ， 而 ] 


有 效 地 保持 了 算法 的 求解 性 


能 09。 

4) 单纯 形 粒 子 群 
optimization, SMPSO) 
在 粒子 群 优化 算法 产生 随机 初始 值 后 ， 使 用 单纯 形 法 对 初 
始 值 进行 处 理 ， 使 用 处 理 后 的 初始 值 进行 粒子 群 优化 算法 的 后 
续 步 骤 。 经 过 单纯 形 法 处 理 后 的 初始 值 如 果 已 经 接近 全 局 最 优 
解 ， 则 经 过 粒子 群 算法 的 迭代 会 很 快 收敛 到 全 局 最 优 解 ， 即 使 
处 理 后 的 初始 值 陷 入 了 局 部 最 优 解 ， 粒 子 群 算法 也 会 有 很 大 的 
概率 使 得 结果 跳出 局 部 最 优 解 。 这 种 处 理 方法 可 以 有 效 发 挥 这 
两 种 算法 的 各 自 优 势 091。 
4.2 数值 实验 数据 

五 种 算法 分 别 对 10 个 30 维 测试 函数 〈 见 附录 ) 独立 进行 
50 次 数值 实验 ， 数值 结果 见 表 1-4， 图 4.2.1-10， 各 算法 的 参数 
如 下 : 
1) DPSO 
POP=3*30+1, ci=c2=1.494, mw 


算法 (19] (simplex method particle Swarm 


a 


设 


=1.4，w,, =0 

2) RPSO 

POP=3*30+1 ,ci=c2=1.494, w=rand(0.5,1) 

3) HSPSO09 

POP=3*30+1 ,ci=0.6, cx=1.6, w=rand(0.5,1), 

变异 系数 lambda=0.85 反射 系数 r=1.75 扩展 系数 e=2.75， 
压缩 系数 c=0.75 

4) SMPSO0I9] 

POP=3*30+1, ci=1.5, c2=2.5, w=rand(0.5,1) 

表 1~3 分 别 记录 了 每 个 算法 50 次 运行 的 最 好 解 、 最 差 解 、 
里 论 准 确 解 之 间 的 误差 ， 表 4 记录 了 每 个 算法 50 次 
的 方差 。 图 6-15 反映 了 50 次 运行 结果 的 平均 值 随 评 
价 次 数 变化 的 情况 .为 了 图 形 对 比 明 显 , 对 纵 坐 标 取 对 数 变换 。 


表 1 最 好 解 与 准确 解 的 误差 

Table 1 Error between the best solution and the exact solution 
函数 名 DPSO HSPSO NM-PSO RPSO SMPSO 
Ackley 1.80E-12 1.82E-07 2.52E-09 1.31E-05 1.8E-06 
Sphere 6.72E-20 1.74E-10 2.36E-13 3.46E-07 1.17E-09 
Rosenbrock 0.591478 7.632082 10.0137 4.090019 4.31444 
Rastrigin 25.86892 14.00768 4.72E-12 18.9097 8.954632 
Griewank 6.51E-14 3.52E-07 5.7E-08 0.000207 1.47E-06 
Styblinski-tang « 28.27347 70.68362 2.89E-05 0.7069 18.24352 
Levy 8.34E-15 4.01E-09 1.09E-11 1.83E-05 1.21E-08 
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Exponential &é 3.33E-16 0 1.11E-16 8.43E-14 5.55E-16 
2 
Zakharov 2.47734 0.00044 0.412443 49.60656 5.132545 i 
Schwefel 2.22 8.42E-05 1.41E-05 3.55E-08 4.03E-05 4.93E-07 100 | 
表 2 最 差 解 与 准确 解 的 误差 i 
Table 2 Error between the worst solution and the exact solution " 
函数 名 DPSO HSPSO NM-PSO RPSO SMPSO 104 上 
Ackley 3.734012 1.98E-06 4.11E-08 0.003062 0.00088 
6 
Sphere 6.95E-13 1.06E-08 1.36E-11 4.3E-05 1.79E-06 > 
Rosenbrock 323.1116 3.03E+03 134.3423 530.7808 118.3099 I 
Rastrigin 94.52081 85.56606 3.054697 77.60666 39.79833 
Griewank 0.08128 0.071211 0.021834 C0.056881 0.061479 图 6 Ackley 函数 平均 值 随 评估 次 数 的 变化 
Styblinski-tang ”226.1875 240.3243 14.13675 169.8959 169.6432 Fig.6 Variations ofmean values with evaluations for Ackley function 
Levy 6.522771 6.546213 1.35E-09 4.365262 0.003875 请 
Exponential 9.33E-15 4.33E-15 4.44E-16 2.65E-11 1.39E-13 
Zakharov 95.01998 126.6204 13.16558 1118.755 96.28544 100 
Schwefel 2.22 0.097455 6.327632 1.51E-07 0.002097 1.14E-05 
表 3 平均 解 与 准确 解 的 误差 a 
Table3 Error between the mean solution and the exact solution 
-10 
函数 名 DPSO HSPSO NM-PSO RPSO SMPSO 
Ackley 0.95083 6.13E-07 1.03E-08 0.000165 5.63E-05 
0 | 2 3 4 5 
Sphere 2.13E-14 3.1E-09 2.02E-12 6.98E-06 2.19E-07 104 
Rosenbrock 57.4121 231.3918 59.97354 64.70934 41.96105 图 7 Sphere 函数 平均 值 随 评估 次 数 的 变化 
Rastrigin 55.61811 44.35636 0.590349 41.96631 17.85171 Fig.7 Variations of mean values with evaluations for Sphere function 
Griewank 0.015881 0.015977 0.000619 0.012978 0.013808 10%0 rosenbrock 2 
Styblinski-tang «138.5399 142.4982 0.346304 90.87972 77.35274 | 
Levy 2.085249 1.308247 2.15E-10 0.413019 7.97E-05 
Exponential &é 1.31E-15 3.82E-16 1.18E-16 2.8E-12 2.7E-14 
Zakharov 17.26885 26.82361 4.285058 402.4986 34.16138 
Schwefel 2.22 0.009718 0.126662 8.22E-08 0.000425 2.73E-06 
表 4 数值 解 的 方差 
Table 4 Variance ofthe solution i 
0 1 2 3 4 5 6 7 8 9 
函数 名 DPSO HSPSO NM-PSO RPSO SMPSO x104 
Ackley 0.808365 1.34E-13 6.14E-17 1.86E-07 2.09E-08 图 8 Rosenbrock 函数 平均 值 随 评估 次 数 的 变化 
Sphere 9.93E-27 5 .46E-18 5.21E-24 8.76E-11 1.81E-13 Fig.8 Variations of mean values with evaluations for Rosenbrock 
Rosenbrock 4878.856 5.21E+05 1000.014 7922.594 796.9568 function 
Rastrigin 203.6541 373.4897 0.651505 173.3405 34.57248 103 ‘en l 
—S— NMPSO 
Griewank 0.00038 0.000301 9.68E-06 0.000178 0.000284 Ss EE 
—A—RPSO 
Styblinski-tang ”1631.403 1426.172 4.162142 1175.893 997.812 村 [ae—oPso | 
SD 0- 和- 哄 6 二 - 太 -6RG 壕 
Levy 2.415501 1.831068 9.04E-20 0.545345 3E-07 "a Bees TT 
Exponential &é 1.96E-30 4E-31 2.22E-33 2.22E-23 1.07E-27 全 AN 
Zakharov 378.6455 1181.804 8.551882 62813.37 480.7299 和 
0 Te 
Schwefel 2.22 0.00041 0.80075 8.06E-16 2.33E-07 4.82E-12 "| Mooeeed 
10 1 2 3 4 5 6 7 
104 


Fig.9 ”Variations of mean values with evaluat 


图 9 Rastrigin 函数 平均 值 随 评估 次 数 的 变化 


ions for Rastrigin function 
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10 ”Griewank 函数 平均 值 随 评估 次 数 的 变化 


Fig.10 Variations of mean values with evaluations for Griewank function 
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图 11 Styblinski-Tang 函数 平均 值 随 评估 次 数 的 变化 
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Variations of mean values with evaluations for Styblinski-Tang 
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12 Levy 函数 平均 值 随 评估 次 数 的 变化 


Fig.12 Variations of mean values with evaluations for Levy function 
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13 ”Exponential 函数 平均 值 随 评估 次 数 的 变化 


Variations of mean values with evaluations for Exponential 
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图 14 Zakharov 函数 平均 值 随 评估 次 数 的 变化 


Fig.14 Variations of mean values with evaluations for Zakharov function 
104 schwefel222 
—S— NM-PSO 
一 HSPSO 
SMPSO 
107 Ta 全 -pso 
一 —— DPSO 
~ 
二 
5 
10 | 网 SL 
Was 
ray re WE 
100 ES al se -ee 4 
ee A 
人 eg 
10 Ws 
| 
1019 
0 1 2 3 4 5 
x104 
图 15 Schwefel 2.22 函数 平均 值 随 评 估 次 数 的 变化 


Fig.15 Variations of mean values with 
function 
4.3 数值 实验 分 析 
4.3.1 表格 数据 分 析 


evaluations for Schwefel 2.22 


观察 表 4，NM-PSO 在 计算 数值 结果 的 方差 上 明显 优 于 其 
他 算法 ， 仅 对 测试 函数 sphere 上 方差 咯 大 于 DPSO， 表 明 NM- 


PSO 在 运行 的 稳定 性 和 结果 的 稳定 性 
察 表 3, NM-PSO 在 运算 结果 的 平均 误差 方 摆 


上 明显 高 于 其 他 算法 。 观 
| 也 处 于 总 体 领先 ， 


对 大 部 分 测试 函数 ，NM-PSO 领先 其 


也 算法 多 个 数量 级 ， 仅 对 


测试 函数 rosenbrock 运算 的 3 


均值 与 其 他 算法 差距 不 大 ， 说 明 
NM-PSO 对 运算 结果 的 平均 误差 优 于 其 他 算法 ， 算 法 的 稳定 性 


强 , 结 果 也 更 精确 .对 表 1 和 2, 对 测试 函数 Rastrigin .Styblinski- 


Tang、Exponential 和 schwefel2.22,， 


其 他 所 有 算法 ,但 对 测试 函数 Rosenbrock, NM-PSO 表现 不 佳 ， 
有 某 些 特殊 的 性 质 ， 观 察 表 2 可 以 得 


这 可 能 是 因为 测试 函数 
出 在 8 个 测试 函数 的 最 差 解 结果 中 ， 


NM-PSO 的 计算 结果 优 于 


本 文 算法 都 处 于 领先 地 位 ， 


特别 在 Rastrigin、Styblinski-Tang、Levy、Exponential .Zakharov、 
Schwefel 2.22 上 领先 了 1 个 数量 级 及 以 上 , 尤 


其 levy 函数 达到 


了 9 个 数量 级 ,体现 NM-PSO 跳出 局 部 最 优 的 能 力 ， 可 以 看 出 


本 算法 的 稳定 性 是 远 远 好 于 划 


他 算法 的 。 对 表格 的 分 析 总 体 表 


明 ，NM-PSO 在 绝 大 多 数 测试 函数 下 性 能 表现 最 佳 ， 即 使 对 表 


现 不 佳 的 测试 函数 ， 性 能 
求解 能 力 总 体 处 于 领先 水 平 。 


不 会 处 于 完全 劣势 ， 所 以 NM-PSO 
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4.3.2 图 表 数 据 分 析 

对 测试 函数 Rastrigin、Styblinski-Tang、Levy 本 文 算法 NM- 
PSO 在 平均 求解 精度 上 相 比 其 他 四 个 算法 表现 具有 压倒 性 的 优 
势 。 并且 NM-PSO 在 跳出 局 部 最 优 解 方 面 展 示 出 了 强大 的 能 力 
在 进行 了 大 约 20000~40000 次 目标 函数 调用 后 ， 
陷入 了 局 部 最 优 解 ,在 之 后 的 调用 中 ,最 优 解 基 本 没有 再 改进 ， 
反观 本 算法 ， 最 终 的 最 优 值 要 好 上 102~104 数量 级 。 对 于 测试 
数 Ackley、Griewank、Exponential、 Zakharov、Schwefel 2.22,， 
本 算法 在 平均 求解 精度 对 其 他 算法 虽然 没有 领先 太 多 ， 但 收敛 
速度 〈 纵 坐标 下 降 速度 ) 全 都 好 于 其 余 四 个 算法 ， 在 同样 的 函 


S| 


数 评估 次 数 下 ， 只 有 HSPSO 算法 在 Griewank 测试 的 前 半 段 性 
能 与 本 算法 不 相 上 下 , 其 余 情 况 NM-PSO 算法 皆 能 够 取得 最 好 


US 


效果 。 对 于 测试 函数 sphere， 三 个 基于 单纯 形 的 粒子 群 算 法 效 
果 均 好 于 RPSO 算法 ， 但 比 DPSO 算法 差 ， 这 可 能 是 单纯 形 与 
粒子 群 算法 融合 结构 决定 的 。 在 对 Rosenbrock 算法 进行 数值 实 
验 时 ， 五 个 算法 均 未 能 很 好 的 达到 全 局 最 优 ， 误 差 都 在 102 数 
量 级 ,更 具体 地 ,对 比 同样 基于 单纯 形 的 算法 HSPSO,SMPSO， 

测试 函数 Ackley 表明 ,三 个 函数 在 最 终结 果 ，, 收敛 速度 方面 均 
要 好 于 传统 粒子 群 算 法 ， 且 其 中 本 文 算 法 依旧 是 最 优 。 对 于 函 
数 Rastrigin、Griewank、Levy、Exponential、Zakharov 三 个 算 
法 中 总 有 两 者 要 好 于 传统 算法 ,体现 了 单纯 形 搜索 带 来 的 好 处 。 


步 比 较 三 者 ， 本 文 算法 NM-PSO 更 是 有 三 者 中 最 佳 的 性 


通过 上 述 对 10 个 测试 函数 的 数值 实验 分 析 得 出 结论 : 在 
分 测试 函数 下 ，NM-PSO 算法 总 体 上 求解 精度 较 高 ， 性 能 
见 优 越 。 


结束 语 


为 了 提高 粒子 群 优化 算法 的 求解 性 能 ， 本 文 一 方面 借鉴 电 
磁 机 制 原理 设计 了 一 种 PSO 迭代 方式 , 另 一 方面 利用 单纯 形 搜 
索 来 加 速算 法 收敛 。 数 值 实验 表明 , 具有 免 参 数 特 征 的 NM-PSO 
能 够 很 好 地 平衡 局 部 搜索 和 全 局 探索 之 间 的 矛盾 。 借 鉴 传统 优 
化 方法 的 成 熟 理 论 、 技 巧 来 指导 设计 进化 算法 是 本 文 的 初衷， 
如 何 设计 一 种 更 好 、 更 简洁 的 算法 融合 方式 是 本 文 今后 研究 的 
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